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■ Abstract. The requirement on the symmetry of the biphoton wave function results 

in a double peak structure of the coincidence spectra of Spontaneous Parametric Down 
Conversion (SPDC) [I]. In this work we report the experimental observation of this 
effect by careful tailoring parameters of the SPDC crystal, the ultra-fast pump laser, 
and the measurement setup. The results are shown to be relevant to quantum state 
' engineering of ultra-fast polarization-, or frequency-entangled mixed states. 
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1. Introduction 

Spontaneous Parametric Down Conversion (SPDC) is a process in which pump photons 
decay in a nonlinear birefringent crystal into pairs of photons of lower frequencies, 
referred to as signal and idler photons (biphotons) (2j, with the phasematching 
conditions being satisfied. The SPDC provides a way for producing entangled states 
of light which are actively used in quantum optics community [3]. The SPDC pumped 
by short laser pulses benefits the timing information of photon pairs and high repetition 
rate which can be useful in practical applications jH|5l|6]. At the same time, in the case of 
short pulses of the pump, the SPDC spectrum becomes wide enough for emitted photons 
to have a wide range of frequencies, which raises the task of measuring idler- and signal- 
photon spectra. A common technique to carry out such investigations is to perform 
a complete set of coincidence measurements with varying frequencies of both emitted 
photons [3 El [9], [10] . The coincidence spectrum of signal photons at a given frequency of 
the idler photons is usually characterized by a rather narrow curve, the width of which 
can be considered as a spectral measure of the correlation of photon pairs [TTJ . However, 
as it was shown in p], in the case of short pump pulses and type-II phase matching under 
appropriate conditions, such coincidence spectrum can be more complicated and can 
have two well separated and well pronounced peaks. This effect can be clearly explained 
by the phase matching conditions and the differences between the refractive indices of the 
ordinary and extraordinary waves in birefringent crystals. However, in theory, the main 
reason of the peak-doubling effect is related to the symmetry requirements imposed 
on the biphoton wave functions. As photons are bosons, and if the SPDC biphoton 
state is pure, then its wave function must be symmetric with respect to the particle 
transpositions \V2\- If the wave function does not appear to be automatically symmetric 
(as in the case of the type-I phase matching [UJ), it has to be completed by the addition 
of terms with transposed particle variables in order to become symmetric. This is the 
way in which a double-peak has been found theoretically jTJ. We confirm the existence 
of the double-peak structure by direct experimental measurements. In accordance with 
the theory, the double-peak structure exists only in the case of pulses shorter than 200 
fs, and under specially arranged conditions of the experiment. 

2. Theory 

Consider the type-II collinear SPDC process pumped by short (femtosecond) laser 
pulses: let both the pump-pulse temporal envelope and its Fourier transform be 
approximated by Gaussian functions. The spectral features of the signal and idler 
photons are determined by their wave function, which will in turn depend on the photon 
frequencies uj\ and U2 as well as the condition that it has to be symmetric with respect 
to the variable transposition 1^2. Such wave function can be presented in the form 

in 

»(.*,<*) « {Q i + Q 2 (l)n^)} , (i) 
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Uq is the central frequency of the pump, r is the pump-pulse duration, and L is the 
crystal length. The two-line columns in Eq. (JTJ correspond to different polarizations of 
photons, (Vj for horizontal and f for vertical polarizations, with the indices 1 and 
2 associated with the numbers of frequency variables u\ and u 2 ; A i2 and A 2 i are the 
phase mismatches 

A 12 (wi, u 2 ) = k p {ui + u 2 , O ) - k Q {ui) - k e (u 2 , O ) 

= n e (uji + uj 2 , O ) — — — - n (wi)— - n e (u 2 , O ) — , (3) 

c c c 

A 2 i(wi, W 2 ) = Ai 2 (w 2 , Ui) = k p (Wi + U) 2 , 0o) - &o(w 2 ) - fc e (Wi, O ) 

= ra e (wi + w 2 , 0o ) n (o; 2 ) n e (cji, O ) — . (4) 

c c c 

Here k ,n and fc e ,n e are wave numbers and refractive indices which correspond to 

ordinary (o) and extraordinary (e) waves in the crystal, and 0o is the angle between the 

crystal optical axis and the direction of photon propagation. The angle 0o is determined 

by the condition of zero mismatch at frequencies where the signal and idler photons 

frequencies are equal to half of the pump central frequency for both terms in Eq. 

([I]): A 12 (w o /2,u;o/2,0o) = A 21 (w /2, u Q /2, O ). Also, A 12 = A 21 at any coinciding 

frequencies U\ = u 2 . But at U\ ^ u 2 mismatches differ from each other (A 12 ^ A 21 ). 

The appearance of two terms in Eq. (TTJ) is related to the obligatory requirement 
that the biphoton wave function must be symmetric with respect to the transposition of 
photon variables 1 #2. An alternative explanation (not contradicting to the first one) 
is based on the observation that in the type-II SPDC process there are two possibilities 
to measure coincidence spectra: the fixed frequency uj 2 can be attributed either to the 
extraordinary (vertically polarized) or to the ordinary (horizontally polarized) waves. 
These two cases exactly correspond to the first and second terms in Eq. ([1]) and to 
two different mismatches A 12 ([3]) and A 21 (jlj). Note that the assumed correspondence 
between polarizations of ordinary, and extraordinary waves in the crystal is determined 
by the conditions of the experiment described in the following section. 

Coincidence spectra of emitted photons are determined by the conditional 
probability densities dw(uj\)/doj\\ W2 , where uj\ is varying and u 2 = const. There are 
two ways of performing coincidence measurements: after being split by a 50/50 non- 
polarizing beamsplitter the down-converted photons are either filtered by polarizers and 
frequency filters in front of the detectors or just by the frequency filters. In the first case 
the spectral conditional probabilities become simultaneously conditional with respect to 
the polarization and are given by squared absolute values of the functions accompanying 
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polarization columns in Eq. ([I]): 
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In the second case, when polarizers are not used, the density matrix from the wave 
function (CD) is constructed and then its trace is taken with respect to the polarization 
variables, which gives 

u 2 ; u)[, u)' 2 ) oc f(ux, u 2 )f(u[, u' 2 ) + f(u 2 , u>i)f(w' 2 , (7) 

where "bar" means averaged over the polarization degree of freedom. Owing to this 
averaging, the density matrix p characterizes a mixed state. Diagonal elements of p 
determine two-frequency probability density 
dw(oJi, uj 2 



p(ui,uj 2 ;uj 1} uj 2 ) oc / (ui,u 2 ) + f (w 2 ,wi), 

QJjJ\QjjJ 2 

as well as the conditional (coincidence) probability density 



(8) 
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As seen from Eq. (J9j), the conditional probability density found for the mixed state (ED) 
coincides with the sum of two probability densities (jSJ) and fl6]) arising in the original 
pure state (0D)-(JI]). There are two ways of obtaining the same polarization-non-selective 
coincidence spectrum. One method consists in using only frequency filters with no 
polarizers. By keeping a given photon frequency u 2 in one channel after the beamsplitter 
and by scanning the filter in the other channel one obtains the spectrum dw^ c >{ujx)/dui 
09]), which is determined by the frequency-dependent mixed state with the density matrix 
p (J7J). An alternative method consists of making two coincidence measurements with the 
polarizer added to the frequency filter in the given-frequency channel, and the polarizer 
oriented vertically in one measurement and horizontally in another one. Each of these 
two measurements will give rise to a single-peak spectrum (JHJ) or (JHJ), but their sum will 
coincide with the double-peak structure determined by Eq. (jUJ). 

FigJTJ shows calculation of the spectral conditional probability density according to 
Eq.()9]) for different pump-pulse durations. All frequency dependencies are recalculated 
to depend on the wavelengths Ao,i,2 = 2^0/^0,1,2 (for convenience of comparisons with 
the experimental data). The parameters of the calculations correspond to the ones 
used in the experiment (see below). For relatively short pump-pulses of r = 110 fs 
(see Fig. la) the spectral conditional probability density has a distinctive "butterfly" 
shape. Such a shape implies that for a fixed wavelength of the signal(idler) photon, 
slightly detuned from the degeneracy, the spectral conditional probability exhibits two 
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Figure 1. Calculation of the spectral conditional probability density for typc-II 
j3— Barium borate (BBO) crystal of the length L = 5 mm, optical axis at 0o = 41.4625°, 
pump central wavelength Ao = 404.7 nm, duration of the pump pulse (a) t = 110 fs, 
(b) t = 250 fs, (c) t = 1000 fs. The energy of the pump pulse is considered to be the 
same for the three cases , i.e. tEq = const, where Eq is the amplitude of the pump 
field. 



distinctive maxima. Each maximum, corresponds to a different polarization component 
of the two-photon state. 

The effect of peak doubling depends crucially on the duration of the pump-pulse. 
The corresponding dependencies of the spectral conditional probability density for the 
case of r = 250 fs, and r = 1000 fs are shown in FigJTb and Fig{Tb, respectively. In the 
case of r = 1000 fs the dependence has only a single maximum at a fixed wavelength 
of the signal (idler) photon, and the effect of peak doubling vanishes. Disappearance of 
the peak doubling in the case of long pulses has a simple qualitative explanation. For 
states with the wave function or density matrix having the form of a product of two 
parts, each depending on two frequencies, widths of conditional spectral distributions 
are determined mainly by the narrower part of the product. In the scheme under 
consideration, these two terms are the pump spectral intensity and the sum of two 
sinc 2 -functions in Eq. (|9]). m the case of short pulses the pump spectrum is wide, 
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and the spectral structure of conditional probabilities is determined mainly by the sum 
of sinc 2 -functions, which gives rise to two peaks. In the opposite case of long pulses, 
the spectrum of the pump becomes very narrow, and the conditional spectra of idler- 
signal photons are determined exclusively by the pump. As the pump spectral intensity 
depends only on the sum of frequencies LO\ + U2 centered at loq, it is clear that at any 
given value of Co> 2 , the conditional spectrum dw^ c \uji) j ' doj\ © has only one peak, and 
the locations of uj\ and U2 are symmetric with respect to w . This symmetry explains 
why the narrow region of large conditional probabilities in FigJIb is concentrated around 
the line tilted at 45° to both the horizontal and vertical axes. 

3. Experiment 

The experimental setup is shown in Figf2J The beam of a mode-locked Ti: Sapphire laser 
(Spectra Physics, Mai-Tai) with a central wavelength of 809.4 nm, pulse duration 80 fs 
at full width half maximum (FWHM), repetition rate 80 MHz, and output power 1.1 W 
is frequency doubled in a 0.3 mm long /3-Barium borate crystal (BBO), placed between 
two confocal lenses with /i=50 mm. The second harmonic (SH) beam is filtered from 
the fundamental beam by low group velocity dispersion dichroic mirrors (not shown), 
and characterized by an autocorrelator (AC, Femtochrome, FR-103XL), and an optical 
spectrum analyzer (OSA, Yokogawa, AQ-6315B). The SH beam has a central wavelength 
404.7 nm, spectral bandwidth 3.5 nm (FWHM), duration 110 fs (FWHM), and power 
200 mW. It is focused by a lens with / 2 =500 mm into a 5 mm long BBO crystal cut for 
type-II (0o ~ 41.45°) collinear frequency-degenerate SPDC at 810 nm. The SH beam is 
vertically polarized (extraordinary wave in the crystal) and its waist inside the crystal is 
200 fim. After the SPDC crystal the SH beam is reflected by a UV-mirror (UVM) whilst 
the SPDC fluently passes through it. In order to compensate for the spatial walk-off in 
the SPDC crystal, a half-wave plate orientated at 45° (HWP1) and two BBO crystals 2 
and 0.5 mm long (total length is equal to the half of the length of the SPDC crystal) with 
their optical axes parallel to the one of the SPDC crystal, are introduced [13J. Then, 
the SPDC passes through another half-wave plate (HWP2) followed by a polarization 
beam splitter (PBS). Such a configuration allows for the routing of different polarization 
components of the SPDC to different output ports of the PBS by rotating the HWP2. 
Next, the SPDC in each output port of the PBS is coupled into a single mode optical 
fiber (SMF) by aspheric lens with ^3=8 mm, placed at the distance 600 mm from the 
SPDC crystal [H]. Each SMF is connected to the input of a home-made spectrometer 
(SI, S2), which consists of an achromatic collimation lens with ^4=60 mm, a diffraction 
grating (Thorlabs, GR-1208) mounted on a motorized rotation stage, and a collection 
lens with /s=175 mm, which focuses the light into a multi-mode optical fiber (MMF) with 
the core diameter of 50 /im. The MMF is connected to an avalanche photodiode (APD, 
Perkin-Elmer, AQRH-14FC) with quantum efficiency 60%. Outputs of the APDs from 
the two spectrometers are addressed to a coincidence circuit (Ortec, Time to Amplitude 
Converter, Model 556) with the time window set to 7 ns. Each spectrometer is tested 
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with a fiber coupled broadband light source, and yields 0.15 nm spectral resolution and 
40% transmission. From the theoretical calculations it is found that slight changes in 
the orientation of the crystal axis (on the order of 0.005°) results in significant changes 
of the observed double-peak structure of the coincidence spectra. To circumvent this 
problem all the presented data has been obtained in a single experimental run without 
realignment of the setup. 
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Figure 2. Experimental setup. The beam of a mode-locked TkSapphire laser is 
frequency doubled, and it is then focused into a type-II BBO crystal, where SPDC 
occurs. The second harmonic (SH) is characterized by an optical spectrum analyzer 
(OSA), and an autocorrclator (AC). UV-Mirror (UVM) filters the SPDC from the 
pump. A half-wave plate (HWP1) and two BBO crystals with parallel optical axes 
(OA) compensate for the spatial walk-off. The SPDC passes through a half-wave plate 
(HWP2), and divided by a polarization beam splitter (PBS). It is coupled into two 
single- mode fibers (SMF), connected to two spectrometers (SI, S2). Each spectrometer 
consists of colhrnating and focusing lenses, and a diffraction grating (DC). Outputs 
of both spectrometers are directed to avalanche photodiodcs (APDs) via multi-mode 
optical fibers (MMFs). APD outputs are addressed to a coincidence circuit. 



4. Results and discussion 

The wavelength (idler), selected by one of the spectrometers (SI), was fixed. In another 
spectrometer (S2) the wavelength (signal) was scanned and the spectral dependence 
of the coincidences was measured for a given polarization component. Then, HWP2, 
initially set at 0°, was rotated by 45° and the above mentioned procedure is carried 
out for another polarization component. The overall coincidence spectrum is obtained 
by summation of the measurements from the two orientations of HWP2. Spectral 
dependencies of the coincidences are tested at four different wavelengths of idler photons 
selected by the spectrometer SI: 809.4 nm, 808.4 nm, 807.3 nm, and 806.3 nm. The 
corresponding results are shown in Fig. |3j 

The experimental data (black squares) demonstrate a fair agreement with the 
theoretical results (red solid lines). With the wavelength of the idler photon detuned 
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Figure 3. Experimental results. Spectral dependencies of the coincidences at different 
fixed wavelengths of an idler photon: (a) 809.4 nm, (b) 808.4 nm (c) 807.3 nm (d) 806.3 
nm. Black points are the experimental data, error bars are the standard deviations, 
red solid lines are the theoretical calculations. 

from the degeneracy, the spectral conditional probability exhibits two distinctive 
maxima. The separation between the two maxima grows with the detuning of the idler 
photon wavelength, in accordance with the theoretically calculated spectral conditional 
probability, using (9). The only free parameter for the fitting of the experimental data 
was the orientation of the optical axis of the crystal. The optimal value was found to 
be 0o = 41.4625°. 

The curve of Fig. la corresponds to the case when the idler wavelength A2 is equal 
to the coinciding central wavelengths of both peaks. Summed together, they produce 
a single-peak curve of a rather unusual shape, with the narrower peak seeming to sit 
on top of the wider one. This combined peak is approximately twice as high as the 
two constituent individual peaks from which this curve is formed. This confirms that 
the observed results correspond to the mixed rather than pure polarization state of 
two photons, because in the case of a pure polarization state, coherent summation of 
individual peaks with coinciding central wavelengths results in a four fold increase of 
the maximal probability compared to that occurring in individual peaks. This is not 
observed. Hence, the experimental results support the theoretical description in terms 
of the mixed polarization state. 
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5. Conclusion 

The presented work shows experimental observation of the earlier predicted effect of 
peak doubling in the coincidence spectra of type-II SPDC and ultra-fast pump pulses. 
The effect occurs only if the pulse duration of the pump is short enough to provide a 
pump spectrum broader than the widest of two peaks. In the case of longer pulses, 
only one of two peaks survives (Fig Jib). The peak-doubling effect is shown to occur 
experimentally only if one analyzes the frequency distributions of photons with their 
polarization degree of freedom being ignored. The experiment has also confirmed that 
the state which is characterized by frequency distributions of photons averaged over 
polarizations is mixed (the peak in Fig. [3b is only two times higher than each of the 
two terms contributing to its formation in Eq. fl9])). The peak-doubling effect is a direct 
consequence of the symmetry requirement of the biphoton wave functions resulting in 
the appearance of two terms in expression (bTJ) for the biphoton polarization-frequency 
wave function. 
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